function [E] = trenergy(stop)
m = 9.109e-28;
E = zeros(stop, 1);
for i=1:1:(stop-1)
    [file,err] = sprintf('tre_dump_%d.dat', i);
    ALL = load(file);
    E(i)= ALL(2)/ALL(1);
end
plot(E,'x');
